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Abstract 

The rare B c — > D s £ + £~decay is investigated by using the most general model 
independent effective Hamiltonian. The general expressions of longitudinal, normal 
and transversal polarization asymmetries for £~ and £ + and the combined asymmetries 
of them are found. The dependencies of the branching ratios and polarizations on the 
new Wilson coefficients are presented. The analysis shows that the branching ratios 
and the lepton polarization asymmetries are very sensitive to the scalar and tensor 
type interactions. These results will be very useful in searching new physics beyond 
the standard model. 
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1 Introduction 



Investigation of rare B meson decays, induced by flavor-changing neutral current (FCNC) 
b —> s,d m transitions, is an important source of new physics. These transitions take place 
at loop level in the SM so they can be used to test the gauge structure of the SM and 
provide a suitable tool of looking for new physics. In rare B meson decays, new physics 
contributions may appear in two different ways; modifying Wilson coefficients in the SM or 
adding new structures in the SM effective Hamiltonian. 

Since the CLEO observation of b — > 57 process [2], decays of B u ^ s mesons have been 
the subject of many investigations. These studies will be even more complete if similar 
studies for B c , discovered by CDF Collaboration [3], are also included. 

In the mean time, the study of the B c meson is by itself quite interesting too, since it 
has some outstanding features [I] [6] . It is the lowest bound state of two heavy quarks (b 
and c) with explicit flavor that can be compared with the charmonium (cc- bound state) 
and bottomium (bb- bound state) which have implicit flavor. The implicit-flavor states 
decay strongly and electromagnetically whereas the B c meson decays weakly. The major 
difference between the weak decay properties of B c and B u ^, s is that those of the latter 
ones are described very well in the framework of the heavy quark limit, which gives some 
relations between the form factors of the physical process. In case of B c meson, the heavy 
flavor and spin symmetries must be reconsidered because both b and c are heavy. 

On the experimental side, like the running B factories in KEK and SLAC, also encour- 
ages the study of the rare B meson decays and most of the rare B c decays are believed to 
be accessible in future experiments at hadronic colliders, such as the LHC-B. The Tevatron 
experiments see around 100 B c semileptonic decays and a luminosity upgrade by a factor 
of ten is discussed for 2013. Atlas and CMS experiments will have more B c but there will 
be more background. This scene may not be hopeful for understanding B c . On the other 
side, rapid progress on experimental techniques are still encouraging. 

Measurement of the lepton polarization is an efficient way in establishing the new physics 
beyond the SM [7] -[To]. In this work we present a study of the branching ratio and lepton 
polarizations in the exclusive B c — > D s £ + £~ [t = fi, r) decay for a general form of the 
effective Hamiltonian including all possible form of interactions in a model independent 
way without forcing concrete values for the Wilson coefficients corresponding to any specific 
model. To make predictions about such an exclusive decay, one requires the additional 
knowledge about form factors, i.e., the matrix elements of the effective Hamiltonian between 
initial and final meson states. This problem, being a part of the nonperturbative sector 
of QCD, lacks a precise solution. In literature there are a number of different approaches 
to calculate the decay form factors of B c — > D s £ + £~ decay such as light front, constituent 
quark models, and a relativistic quark model proposed in |16j . In this work we will use the 
weak decay form factors calculated in [16]. 

The work is organized as follows. In section 2, we derive the matrix element starting from 
the effective Hamiltonian for the quark level process and using the appropriate form factors. 
Then, we present the model independent expressions for the longitudinal, transversal and 
normal polarizations of leptons and combined lepton-antilepton asymmetries. We give our 
numerical results and discussion in section 3. 
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2 Effective Hamiltonian and Lepton Polarizations 



The B c — > D s £ + £~ decay is described at the quark level by the b — > s£ + £~ transition in 
the standard effective Hamiltonian approach. This Hamiltonian includes all possible terms 
calculated independent of any models. The effective Hamiltonian for this process can be 
written in terms of twelve model independent four- Fermi interactions, as follows [T7j: 

H eff = ^V ts V t t\C SL sia^Lbl^£ + C BR sia^Rb£Yi 

V2vr [ q q 

+ cti sli^l i Ll n L + cl% s LllM b L l Rl n R + c RL s Rlll b R l Ll n L 

+ C RR sbJ^r In^in + C LRLR s L b R £ L £ R + C RLLR s R b L £ L £ R (1) 
+ C LRRL s L b R £ R £ L + C RLRL s R b L £ R £ L + C T sa^b £a^£ 

+ iC TE e^sa^bIa a/3 £^ , 

where L = 1 — 7 5 /2 and R = 1 + 75/2 are the chiral projection operators and Cx are the 
coefficients of the four-Fermi interactions. The coefficients Csl and C BR are the nonlocal 
Fermi interactions corresponding to —2m s C^^ and —2mbCj^ in the SM, respectively. The 
Cll, Clr, C R l and C RR terms are the vector type interactions, two of which are vector 
interactions containing C l £ L and C^ R do already exist in the SM in combinations of the 
form (C 9 e// - C 10 ) and (C 9 e// + C l0 ). Therefore, we write 

^LL — °9 ^10 1 ^LL 1 

ritot sieff \ r 1 _l_ C 

^LR — °9 ' t-^lO T t-'LR ; 

so that that C l fl and C l f R describe the sum of the contributions from SM and new physics. 
The terms with coefficients Clulr, C R ll R , Cl RR l and C R l R l describe the scalar type 
interactions and the last two terms, CV and Cte, describe the tensor type interactions. 

Having the general form of four-Fermi interaction for the b — > s£ + £~ transition, the next 
task is to calculate the matrix element for the B c — > D s £ + £~decaj which can be expressed 
in terms of the invariant form factors over B c and D s . These form factors are weak decay 
form factors [16] • 

(D s (p Ds ) \sia^q v b\ B c (p Bc )) = ^ \(P Bc + P Ds ) q 2 - q^m 2 Bc - m 2 Ds )] , (2) 

m Bc + m Ds 1 r c j 

(D s (p Ds ) \s^b\ B s (p Bc )) = F + {q 2 ){P Bc + P Ds )^ + F^q 2 )q, , (3) 

(D s (p Ds )\sb\B s (p B J) = F.—^ + F ^'^' (4) 

m b — m s mb — m s 



and 



(D s (p Ds ) \sa^b\ B s (p Bc )) = 1 \{P Bc + P Ds )^q u - q,{P Bc + P D J 1 , (5) 



m Bc + m D 

where q = p Bc — po s is the momentum transfer. 
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We can write the matrix element of the B c — > D s £ + £ decay using Eq. ([I])-© as 

Ga 



M(B C 



D a 



V tb V*{£Yi A(Pi 



p, 



D s )n 



Bq, 



+£lV [C(P Bc + PdX + + U Q + hstN 
+Ua^£(iG) [(P Bc + P Ds ),q v - (P Bc + Pd s )u + q, 

+Ua ol pie^H\{P Bc + PdXQu - (Pb c + P Ds )»q» 



(6) 



where 



.4 

B 

C 
D 

Q 

N 
G 
H 



{CTl + C? R + C RL + C RR )F + (q 2 ) - 2(C SL ) + C BR ) 



(Ci°L + Cf L + C RL + C RR )F4q 2 ) + 2(C S L + C BR )F T 

(Clr + Crr - Cll - C RL )F+(q 2 ) , 
{Clr + Crr ~ Cll ~ C RL )F^(q 2 ) , 



m Bc + m Ds 
(m Bc - m Ds 



q 



LRLR 



+ c 



RLLR 



+ c. 



LRRL 



+ c 



RLRL ) 



q 



m Bc — m 



[ClRLR + C R LLR — Clrrl — Crlrl) 

Prp 



"Mb — m s ™b — m s 

9 2 ^ 

F q +F + ^ 



m 



m b - m s 



-Cn 



-G 



m Bc + m Ds 



TE 



m Bc + m Ds 



(7) 



In order to calculate the final lepton polarizations, we define the orthogonal unit vector 
S~[" in the rest frame of i~ and <S$ in the rest frame of £ + and the polarization of the 
leptons along the longitudinal (L), transversal (T) and normal (N) directions, as done 
before [ZKHj, by 



Q-fJ, 



Q-fJ, 



St." 



q+rl 

°N 



J7 1 



(o,e L -; 

(o,ef ; 



(0,e L + ) 



(0, e+) 




0. 



V- 



' \p-\J 

p x p_ 
\p x p_ 

0, e N x e L 
P+ \ 



0. 





(0,e 



\p+\J 
p x p + 

' \p x P+\ 



0, ejy x 



(8) 



3 



where p± and p are the three momenta of £ and D s meson in the center of mass (CM) 
frame of the lepton pair system, respectively. The longitudinal unit vectors S£ and S£ are 
boosted to the CM frame of £ + £~ by Lorentz transformation, 



L, CM 



C+IJ- 
°L,CM 



\P-\ EeP- 



m e m e \p. 



E t p. 



m e 



me \p. 



(9) 



while vectors of perpendicular directions are not changed by boost. 

Asn ± being any spin direction of the £ ± , in the rest frame of the leptons, the differential 
decay rate of the B c — > D s £ + £~decay can be written in the following form: 



dT(n d 
ds 



1 (dT s 

2 \~ds~ 



1 + I P±e£ + P^e± + P±e 



N^N 



rp K^rp 



n 



(10) 



Here, s = q 2 /m 2 Bc , the superscripts + and correspond to the £ + and £ cases and (dT/ds)o 
corresponds to the unpolarized decay rate, whose explicit form is 



'dT 
ds 



G 2 a 2 m Bc 

2 14 7T 5 



V tb V*\ 2 V\vA 



(11) 



where 



A 



1024 



sv 2 \m% \H\ 2 + 



256 



s(3-2v 2 )\m Bc \G\--(r 



(v 2 -3)\m B \A 



l28me\m Bc Re(AG*) + 32(1 
, 2 C 4 



r m, 



'm 2 B Re{CD*) + 16sm e m 2 Bc \D\ 



+ Asm 2 Bc |Af + jm Bc s[2X - (1 - v 2 )(2\ - 3(1 - r) 2 )] \C\ 2 
+ 16(1 - r)m e m% Re(CN*) + 16sm^m| Re(DN*) + 4sv 2 m 2 B \Q\ : 



(12) 



and A = 1 + r 2 + s 2 — 2r — 2s — 2rs, r = m 2 D /m 2 B and lepton velocity is v 
The polarizations Pf;, P? and P N in Eq. (TlOl) are defined by 



1 — 4m 2 / q 2 



dT 

dq 2 



(n 



dT 

dq 2 



[n 



—e- 



dT 

dq 2 



[n 



± 



-e- 



for i = L, N, T. Here, P^ and P^ represent the longitudinal and transversal asymmetries 
of the charged lepton £ ± in the decay plane and P N is the normal component to both of 
them. After calculations, the longitudinal polarization of the £ is 



pi 



\m\v 



A 



2 1 6 

=F -m 2 B \Re{AC*) ± —m 2 B m e XRe{CG* 

3 3 

128 



32 



m 2 B me\Re{AE" 



-2me(l - r)Re(CQ*) + —m B \sRe{GH*) - 2m e sRe(DQ*) - sRe(NQ* 



(13) 



where A is given in Eq. (fT2l . 

In a similar way, we find the transverse polarization 



p± _ B c V 



T 



2 32 
± -m/l - r)ite(AC*) ± — m?(l - r)Re(CG*) =f 2mgRe(AD* 
s s 



±32m 2 e Re(DG*) T Re(AN*) ± 16m e Re(GN*) + v 2 Re(CQ*)\ , (14) 

In the limit of me — > 0, the transverse polarization is due the scalar terms. This can give 
new information about new physics. 

Finally, the normal polarization P N is given by 

L! ' ' " 2m e Im{CD*) - Im(CN*) ± Im(AQ*) =f 16m e Im(GQ*) . (15) 



p± 

A 

In this work we assume all form factors and all new Wilson coefficients are real. Therefore, 
the only contribution to Pfj in Eq. (|T5|) comes from ±/m(AQ) term since only the function 
A has an imaginary part coming from Cg . For this reason, P^ = P^ = in the SM 
and scalar terms in Q makes normal polarization nonzero beyond the SM. This observable 
result gives useful clue about new physics. 

Since in the SM P L + P£ = 0, P T - P^ » and P N + P+ = 0, combined analysis of 
the lepton and antilepton polarizations can be another useful source of new physics [T7] . 

Using Eq. ( fT3l) we get combined longitudinal polarization 

" n 128sAm| Re{GH*) - 32m t m\ XRe{AH*) 



p— I p+ Be 

Fl +n ~^A~ - 

-3sRe(NQ*) - 6(1 - r)m t Re(CQ*) - 6m^i2e(.DQ*)] , (16) 

and combined transversal polarization is the difference of the lepton and antilepton polar- 
izations and can be calculated from Eq. ffT4l) 

p-_ P + 4m| c 7rVsA r 32m|(l - r) n , 2m/ (1 - r) 



T ±T A 



-Re(CG*) + — ^ '-Rt{AC 



s s 

,2; 



+2m £j Re(AD*) - 32mjRe(DG*) + Re(AN*) - 16m £ Re(GN*)\ . (17) 

The terms containing the SM contribution to the P£ + P£ in Eq. (lT6"j) . completely cancel 
so that any nonzero measurement of this value in the experiments will provide essential 
evidence of new physics beyond SM. The combined normal polarization, P^ + P^, is zero 
beyond the SM, since P^ and P£ receive contribution only from Im(AQ) term with opposite 
sign. 

A last note before going into details of numerical analysis is in order. In the expressions 
of the lepton polarizations we note that they all depend on s and the new Wilson coef- 
ficients. Because of experimental difficulties of studying the polarizations of each lepton 
depending on both quantities, it would be better to eliminate the dependence of the lepton 
polarizations on s, by considering the averaged forms over the allowed kinematical region. 
The averaged lepton polarizations are defined by 



[Pi) 



(l-m D Jm Bc ) 2 d13 

Pi—ds 

(2 mt/m Bc ) 2 as 

(l~ mD Jm Bc ) 2 

—ds 

(2m e /m Bc ) 2 ds 
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3 Numerical analysis and discussion 



Before going on our numerical analysis of the branching ratios and the averaged polarization 
asymmetries < P£ >, < P^ > and < P^ > of i~ for the B c — > D s £ + £~ decays with i = fi,T 
as well as the lepton-antilepton combined asymmetries < P£ + P£ > and < P^ — P^ >, 
let us first introduce the input parameters used in this work: 

m Bc = 6.50 GeV, m Ds = 1.968 GeV, m b = 4.8 GeV, m M = 0.105 GeV, m T = 1.77 GeV, 
\V tb V*\ = 0.0385, a~ l = 129, G F = 1.17 x 10" 5 Ge1/" 2 , 

t Bc = 0.46 x 10~ 12 s, Cj ff = -0.313, C 9 e// = 4.344 , C 10 = -4.624 (19) 

Details of the values of the individual Wilson coefficients in the SM at fi ~ mj scale can 
be found in [15] . 

The given value of corresponds only to the short-distance contributions, but we 
know that Cg also receives long-distance contributions due to conversion of the real cc 
into lepton pair £ + £~ , and they are usually absorbed into a redefinition of the short-distance 
Wilson coefficients: 



c 9 e// (/i) = c 9 (/i) + r(/i) , (20) 



where 



Yfa) = Y reson + h(y,s)[3C 1 (fi)+C 2 (fi)+3C 3 (fi) + C 4 ( f i) + 3C 5 (fi) + CM] 
s) {ACM + 4C 4 M + 3C 5 (/x) + CM) 
- ±h{0,8)[CM+3CM) (21) 
+ ^(3C 3 (/i) + C4(/i) + 3C 5 (/i) + C 6 (/i)) , 

and ?/ = m c jm b , and the functions s) arises from the one loop contributions of the four 
quark operators 0\,...,Oq explicit forms of which can be found in |18j - [20] . Parametrization 
of the resonance cc contribution, Y reson (s), given in Eq. (l2TI) can be done by using a Breit- 
Wigner shape with normalizations fixed by data given in [21] 

_ 3 v TrTQj -> l+r )m v< 

-< resort^ 9 7, o 1 • ri 

a em v % t sm 2 B - m Vt + im Vi T Vi 
x [(3C 1 (fx) + C 2 ^) + 3C 3 (fx) + CM + ^CM + CM)}, (22) 

where the phenomenological parameter k is usually taken as ~ 2.3. 

The new Wilson coefficients are the free parameters in this work, but it is possible to 
establish ranges out of experimentally measured branching ratios of the semileptonic rare 
B-meson decays 

BR(B -> Ktt) = (4.8±J;g ± 0.3) x 10~ 7 , 
BR(B -> K* n + yT) = (1.27+°;^ ± 0) x 10" 6 , 
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reported by Belle and Babar collaborations [22] 
rare B-decays in the B° — > fi + fi~ mode [23]: 



[23] . also the upper bound of pure leptonic 



BR(B G 



fi fi 



< 1.5 x 10 
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All new Wilson coefficients are taken as real and varying in the region —4 < Cx < 4, 
compliant to this upper limit and the above mentioned measurements of the branching 
ratios for the semileptonic rare B-decays. 

The new Wilson coefficients in Eq.([T]), the helicity-ffipped counter-parts of the SM oper- 
ators, Crl and Crr, vanish in all models with minimal flavor violation in the limit m s — > 0. 
However, in some MSSM scenarios there exist finite contributions from these vector oper- 
ators even for a vanishing s-quark mass. In addition, scalar type interactions can also 
contribute through the neutral Higgs diagrams, multi-Higgs doublet models and MSSM, 
for some regions of the parameter spaces of the related models. In literature there are some 
studies to establish ranges out of constraints under various precision measurements for these 
coefficients (see e.g. [25]) and our choice for the range of the new Wilson coefficients are in 
agreement with these calculations. 

To make some numerical predictions, we also need the explicit forms of the form factors 
F + , F_ and Ft- Since there are two heavy quarks the non-relativistic effects and also the 
additional contributions from hard interactions might give useful information. In this work 
we have not considered these effects and used the results of [16], calculated in a relativistic 
constituent quark model in which q 2 dependencies of the form factors are given as 



F(q 2 



F(0) 



2 ! 



(l — as + bs 2 ^ 

where values of parameters F(0), clf and &f for the B c — > D s decay are listed in Table 1. 





F(0) 


a 


b 


F + 


0.186 


2.48 


1.62 


F_ 


-0.190 


2.44 


1.54 


Ft 


0.275 


2.40 


1.49 



Table 1: B c meson decay form factors in a relativistic constituent quark model without 
impulse approximation . 

Before the discussion of the results of our analysis given in a series of figures, we give our 
SM predictions for the longitudinal, transverse and the normal components of the lepton 
polarizations for B c — > D s £ + £~ decay for fi (r) channel for reference: 



<Pl> 
<P T > 

<P N > 



-0.8457 
-0.0948 
0(0). 



-0.1774) , 
-0.6189) , 
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In Figs. ([T]) and (J2]), we give the dependence of the integrated branching ratio (BR) on 
the new Wilson coefficients for the B c — > D s fi + yU~and B c — > D s t + t~ decays, respectively. 
From these figures we see that BR depends strongly on the scalar and tensor interactions 
and weakly on the vector interactions. It is also clear from these figures that dependence 
of the BR on the new Wilson coefficients is symmetric with respect to the zero point for 
the muon final state, but such a symmetry is not observed for the tau final state for the 
tensor interactions. Another remark is that, for muon case Cte is dominant while for tau 
case Ct becomes more dominant. 

In Figs. ([3]) and (HI), we present the dependence of averaged longitudinal polarization 

< P£ > of £~ and the combined averaged < P£ + P£ > for B c — > D s ji + pT decay on the new 
Wilson coefficients. It is observed that the dominant contribution for < P£ > comes from 
the scalar interactions of the type Clrrl and Crlrl which are identical and symmetric with 
respect to Cx = while the combined averaged < P£ + P£ > is sensitive to that of scalar 
type interactions only. It is a well-known fact that vector type interactions are canceled 
when the longitudinal polarization asymmetry of the lepton and antilepton is considered 
together. This is the reason < P£ + P£ > does not exhibit any vector type dependence in 
Fig. (j3J). It is also interesting to note that < P£ + P£ > is positive for Clrrl and C crlrl 
and negative for other scalar contributions. In addition they are symmetric with respect to 
C x = 0. 

Figures (jHJ) and (jEJ) are the same as Figs. ([3]) and (BJ, but for B c — ■> D s t + t~ . < P£ > 
strongly depends on scalar type interactions and also sensitive to tensor type interactions. 
In case of vector interactions, they are nearly identical for Cx > and keeping their SM 
values. As in the muon case, < P£ + P£ > depends on scalar interactions and Cte tensor 
interactions. Cte effect is comparably higher than that of muon case and it is negative 
(positive) for Cx < (Cx > 0). For < P£ + P£ > in both \i and r decays, there is no SM 
effect and any nonzero experimental results of < P£ + P£ > will be the evidence of new 
physics. 

In Figs. (Ej) and (jSJ), we present the dependence of averaged transverse polarization 

< > of i~ and the combined averaged < P^ — Pj. > for B c — ► D s £ + £~ decay on the new 
Wilson coefficients. The significant effect of scalar Clrrl and Crlrl can be seen from Fig. 
([7]). Here, < > is negative (positive) for Cx ~ (Cx ~ 0). As a last remark, all other 
contributions are negative except Clrlr and Crllr scalar terms for Cx ~ 0.8. Unlike the 

< P T >, Clrlr and Crllr make greater contribution to < P T — P^ > as seen in Fig. 
(JHJ). Considering Figs. (JTj) and © together, determination of the sign and magnitude of 
the scalar observables can also give useful information about the existence of new physics. 

Figures and flTDJ) are the same as Figs. ([7D and but for B c — > D s t + t~ . In both 
figures, the effects of scalar contributions are clear. Considering < P^ — PjT > it can be 
noted that values of obtained data are two times that of < P^ > values. 

Figures ffTTl) and ffl2|) . we present the dependence of averaged normal polarization < 
Px > of l~ for B c — ► D s fi + fi~a.nd B c — > P s r + r _ decays on the new Wilson coefficients. 
Normal components of polarization of both channel depend only on scalar type interactions. 
As discussed before, SM contribution of normal polarization is zero so any measurable values 
should come from new physics. These values will also give promising information on sign 
and magnitude of < P^ > because < P^ > is negative (positive) for Cx < 0(Cx > 0). 

In conclusion, using the general, model independent form of the effective Hamiltonian, 
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we present the most general analysis of the lepton polarization asymmetries in the rare 
B c — > D s £ + £~ decay. The dependence of the longitudinal, transversal and normal polar- 
ization asymmetries of £~ and their combined asymmetries on the new Wilson coefficients 
are studied. The lepton polarization asymmetries are very sensitive to the existence of the 
scalar type interactions and in some cases tensor type interactions worth to be considered. 
Individually, Clrrl and Crlrl play a significant role throughout this work. In all types of 
analysis the following terms are found identical: Crr = Clr, Crl = Cll, Clrrl = C^lrl 
and Clrlr = C^llr- Moreover, in the most cases polarization effects change their signs as 
the new Wilson coefficients vary in the region of interest, which is useful to determine the 
sign in addition magnitude of new physics effect. A last note on combined asymmetries, a 
well known SM result that < P L + P£ >= 0, < P^ + P^ >= and < P T - P^ >~ 
in the limit rri£ — > 0. Therefore any deviation from these relations and determination of 
the sign of polarization is decisive and effective tool in searching new physics beyond the SM. 
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Figure 1: The dependence of the integrated branching ratio for the B c — > D s fi + {i~ decay 
on the new Wilson coefficients. 
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Figure 4: The dependence of the combined averaged longitudinal lepton polarization 
< P£ + P£ > for the B c — > D s fi + fi~ decay on the new Wilson coefficients. 
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Figure 5: The same as Fig. but for the B c — > D s t + t decay. 
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Figure 10: The same as Fig. (JED, but for the B c — > D s t + t decay. 
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Figure 11: The dependence of the averaged normal polarization < P N > of £ for the 
B c — > D s ji + fi~ decay on the new Wilson coefficients. 
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